The static potential in Jv =4 supersymmetric Yang-Mills at weak coupling 
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We compute the static potential associated to the locally 1/2 BPS Wilson loop in M=A super- 
symmetric Yang-Mills theory with 0(X 2 /r) accuracy. We also resum the leading logarithms, of 
0(\ n+1 In" A/r), and show the structure of the renormalization group equation at next-to-leading 
order in the multipole expansion. In order to obtain these results it is crucial the use of an effective 
theory for the ultrasoft degrees of freedom. We develop this theory up to next-to-leading order 
in the multipole expansion. Using the same formalism we also compute the leading logarithms, of 
C(A 71+3 In" A/r), of the static potential associated to an ordinary Wilson loop in the same theory. 
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There is a huge interest in the study of the Yang- 
Mills theory with J\f = 4 supersymmetry in four di- 
mensions. One of the reasons is the conjectured ex- 
istence of a correspondence between the 7V=4 super- 
symmetric Yang-Mills and the type IIA superstring 
theory on an AdSs x S 5 background This du- 
ality is known as the AdS / CFT correspondence and 
has importance consequences, since it allows to com- 
pute the strong 't Hooft coupling limit (A = N ? g ) 
of certain correlators in the A/"=4 supersymmetric 
Yang-Mills theory for large N c . This is so because 
this limit becomes equivalent to the classical limit 
on the string theory side for which computational 
techniques exist. 

Checks of this conjecture are difficult to obtain, 
since the only quantitative approach to gauge theo- 
ries is based on computations at weak coupling. In 
some cases it is possible to check the conjecture with 
the result at weak coupling. This usually happens 
when non-renormalization theorems exists that per- 
mit to perform the computation on the perturbative 
side exactly. 

In other cases no such checks exist and, usually, 
one can only study both the weak and strong cou- 
pling limit with increasing degree of accuracy hoping 
to gain further input on how the extrapolation from 
weak to strong coupling takes place. In this paper 
we concentrate in one of such correlators, the locally 
1/2 BPS static Wilson loop, or, more specifically, the 
large time limit of its logarithm. Its strong 't Hooft 
coupling limit for large iV c has been computed us- 
ing the AdS/CFT correspondence in 0, H|. O n the 
other hand the question whether the understanding 
of the infrared structure of the static potential in 
the weak coupling regime may shed some light on 
the AdS/CFT correspondence has been addressed 
in Refs. [1, HI]. In these references some infrared 
divergences were found, which allowed the authors 
to obtain the leading logarithmic correction to the 
tree level result. These infrared divergences have 
a similar origin to those found in the QCD static 



potential at weak coupling [6], which, however, in 
this case first appear at three loops. They are due 
to the existence of degrees of freedom with energy 
of C(A/r), which we will call ultrasoft in what fol- 
lows. This scale is much smaller than the soft scale 
~ 1/r at weak coupling. In any case, it is some- 
what surprising that, after so much work, the static 
potential has not even been computed with next- 
to-leading-order (NLO), ie. 0(X 2 /r), accuracy yet. 
The reason is that, at this order, ultrasoft effects 
enter into the game and an infinite resummation of 
diagrams is needed in order to obtain the desired 
accuracy. However, this problem can be revisited 
from an effective field theory perspective, as it has 
already been done for the QCD case 0. As we will 
see, by doing so, the problem trivializes and we will 
able to compute the static potential with C(A 2 /r) 
accuracy. The use of effective theories will also al- 
low us to write renormalization group equations for 
the static potential. By solving them we will also 
obtain the static potential with leading- log (LL), ie. 
C(A™ +1 ln"A/r), accuracy. 

We will also use this example to give full details of 
how the factorization between the soft and ultrasoft 
scale takes place in the static potential for a specific 
computation including finite pieces. In the case of 
QCD this factorization takes places at three loop 
and the full computation does not exist yet. 

Finally, using the same formalism we will also 
compute the static potential associated to an ordi- 
nary Wilson loop with NNLL accuracy up to the 
two-loop matching condition. 



I. THE STATIC SINGLET ENERGY 



The Af = 4 SUSY Lagrangian reads 
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<I> and \& represent a scalar and majorana parti- 
cle respectively. The adjoint covariant derivative 
reads: D^-f = <9 M (-) a - g f abc 'A^-f , i range from 
1 to 6, and the dots represent interactions be- 
tween the scalars and the majorana particles or self- 
interactions between the scalars. 

We now take the locally 1 /2 BPS Wilson loop (for 
a motivation of this definition see Q) 



Wr 



N, 
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where $„ = $ • n, n is a six-dimensional vector with 
n 2 = 1, and C represents the path followed by the 
Wilson loop. In our case we consider a static Wilson 
loop and its associated singlet static energy 



E 8 (r)= lim -ln(Wb) 
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(3) 



where Wo is the rectangular Wilson loop with 
edges xi = (T/2,r/2), x 2 = (T/2,-r/2), Vl = 
(-T/2,r/2) and y 2 = (-T/2, -r/2). The symbol 
( ) means the average over the massless fields. 

In order to describe the static Wilson loop at the 
dynamical level we consider the Lagrangian ([1]) phis 
static sources in the appropriate representation 



(4) 



where 



C 



= ^(ido-gA -g^ n )ij+xi(ido+gA^-g^)xc 

ip and Xc (the conjugated field of \) correspond 
to static sources in the fundamental and anti- 
fundamental representation respectively. The case 
with only one static source has been considered in 
Ref. [1]. 

With this Lagrangian the static singlet energy can 
be obtained from the following Green function 



(0\ X \x 2 )W(x 2 ,x 1 )^{x 1 ) 

x^(y 1 )W(y 1 ,y 2 )x(y2)\0). 
S a ( Xl - yi )5 3 ( X2 -y 2 )(W a ), 



(6) 



where W(x 2l x\) = Wc for a straight path C with 
initial and final points Xi and x 2 respectively. 



shown to be so for the QCD static potential 0, El • 
In that case the effective theory was pNRQCD [ijj. 
For a review see [HJ ■ Here we would like to construct 
the Af = 4 supersymmetric version of pNRQCD in 
the static limit. We will do so at next-to- leading or- 
der in the multipole expansion. As we will see the 
main difference is the existence of massless scalars. 
On the other hand the heavy quark and antiquark 
rearrange in a singlet or octet configuration under 
(ultrasoft) gauge transformations as in QCD. The 
simplification arises from the existence of two dis- 
parate scales: the soft ~ 1/r and the ultrasoft ~ A/r 
scales, and that we only aim for a description of the 
dynamics at energies of order ~ A/r. Therefore, de- 
grees of freedom with energy ~ 1/r can be integrated 
out and one can factorize the physics associated to 
each scale. Note that we live in the opposite limit 
to strong coupling. For us E s (r) <C 1/r, whereas in 
the strong coupling limit E s (r) 3> 1/r. 

Integrating out the soft scale, 1/r, from ((J) we are 
left with an effective theory where only ultrasoft de- 
grees of freedom remain dynamical. The surviving 
fields are the ip~x states (with ultrasoft energy) and 
the massless ultrasoft gluons, scalars and fermions. 
The ijj-x states can be decomposed into a singlet 
(S) and an octet (O) under colour transformation. 
The relative coordinate r = xi — x 2l whose typical 
size is the inverse of the soft scale, is explicit and 
can be considered as small with respect to the re- 
maining (ultrasoft) dynamical lengths in the system. 
Hence the massless fields can be systematically ex- 
panded in r (multipole expansion), and the effective 
Lagrangian is constructed order by order in r. As 
a typical feature of an effective theory, all the non- 
analytic behaviour in r is encoded in the matching 
coefficients, which can be interpreted as potential- 
like terms. 

In order to have the proper free-field normaliza- 
tion in the colour space we define 



S = 



/Nr. 
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where Tp — 1/2. 

The effective Lagrangian density that can be con- 
structed with these fields and that is compatible with 
the symmetries of Eq. @ is given at the next-to- 
leading order in the multipole expansion by: 

£ us = TrJst (id -V s (r) + ...)S 



II. EFFECTIVE THEORY 

The use of effective field theories allows us for a 
efficient description of the infrared structure of the 
static potential at weak coupling. This has been 
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+gV A (r)Tr {C^r -ES + S f r-EO} 
+.9^Tr{0t{r-E,0}} 

_ 5 Ml T r{Ot[r.(DO>„),0]}, (8) 

where R = (xi + x 2 )/2, S = S(r,R,t) and O = 
0(r,R, t) are the singlet and octet wave functions 
respectively. All the gluon and scalar fields, as well 
as the derivative of them, in Eq. §8§ are evaluated 
in R and t. In particular F*° = E = E(R, t) and 
iD O = id a O-g[A (R,t), O]. V x (r) are the match- 
ing coefficients of the effective Lagrangian. They are 
determined by matching the effective and the un- 
derlying theory at a scale v smaller than 1/r and 
larger than the ultrasoft scales. Since we are at 
weak coupling, the matching can be done pertur- 
batively. At the lowest order in the coupling con- 

2 

stant we get a>v B = «v = ot s = £^> Va = Vb = 
V$ = V$ = I, where we have defined (Ca = N c . 
C F = (A 2 - l)/(2iV c )) 

V s (r) = -2C F ^, (9) 

which correspond to the singlet and octet heavy ip~x 
static potential respectively. 

Note that we distinguish between the singlet static 
energy E s (r) and the singlet static potential V s (r). 
The reason for that has been discussed in detail in 
Ref. 0. V s (r) corresponds to the matching coef- 
ficient that appear in the effective theory and it is 
the proper object to be introduced in a Schroedinger 
equation in case we were working with particles with 
large but finite mass. On the other hand E s (r) di- 
rectly corresponds to Eq. ([3]) and represents the 
energy of two static particles (one could do a similar 
distinction for the octet static potential and energy, 
which in the QCD case corresponds to the hybrid 
energy). 

Charge conjugation proves to be a very useful 
symmetry to eliminate operators in the effective La- 
grangian. The effective Lagrangian is invariant un- 
der charge conjugation plus particle <-> antiparticle 
exchange. In particular singlet, octet, gluon and 
scalar fields transform as: 

S(*,r,R) ^<T 2 S(f,-r,R) T cr 2 , (10) 

0(i, r, R) -> a 2 0{t, -r, R) T cr 2 , 

A^t, R) -> -A^t, R) T , $ n (t, R) -> $ n (t, R) T . 

Asking for invariance under the above transforma- 
tions the following operators are eliminated: 

< 5£ = -. 9 Vz 3 Tr{O t [$„,0]} , (11) 



SC = -gV E r (Tr {S t (D$„)0} + h.c.) , (12) 

5C = - 9 V F rTr {0* {(D$„), O}} . (13) 

One can also easily check perturbatively up to order 
a 2 that the matching coefficient of these operators 
are zero. Therefore, they do not affect the results 
obtained in this paper for the singlet static potential 
and energy. 

We have to mention that we have not explored 
all the constraints supersymmetry (or any other 
remaining underlying symmetry) may impose on 
the matching coefficients of the effective Lagrangian 
and, accordingly, on the structure of the renormal- 
ization group equation. This would need a dedicated 
study that goes beyond the aim of this work. In the 
case of QCD some analysis have been performed for 
the underlying Poincare symmetry of pNRQCD in 
Ref. [3. 

Finally, let us note that the octet singlet potential 
is 1/iV 2 suppressed compared with the singlet poten- 
tial. Therefore, it can be neglected in the large N c 
limit. 



III. RG 

The renormalization group equation of the singlet 
static potential has the following structure 

v^-V 8 = 7i(K - V ) + - fy r 2 (V s -V Q f + ---, (14) 
av 

where v is the factorization scale, and the dots refer 
to higher orders in the multipole expansion. This 
structure follows from the multipole expansion and 
the mass gap between the octet and singlet static po- 
tential. The anomalous dimensions 7* have structure 
themselves. 71 can be computed as an expansion in 
a s (note that some powers of ^Jal get associated a 
power of Vq> or V# ). For 73 we also have depen- 
dence on Vjl (but in a very specific way) : 

73 = V%(2C F a a )F(a,,V*,V* ) (15) 

For both anomalous dimensions we can easily obtain 
the lowest non-trivial contribution. They read 

2 2Cf£^ +0(q5)> (16) 

7T 

l. = \ 2 -^Vl+0{al). (17) 

They come from the computation of the ultraviolet 
behavior of the diagram in Fig. [I] with the V$ or V A 
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vertex. Note that the computation of 71 is explic- 
itly gauge invariant in the effective theory and only 
needs the computation of one diagram. Therefore, it 
does not need of the delicate cancellation between in- 
frared divergences that takes place at the soft scale, 
in particular, in the evaluation of the leading loga- 
rithm made in Ref. upon which the contribu- 
tion of each particular diagram is gauge dependent. 
Quite remarkably, the leading order contribution to 
73 is the same as the one one would obtain in pure 
QCD, since it is only due to the chromoelectric field. 
On the other hand in QCD 71 = . 



FIG. 1: One loop contribution to the singlet propagator. 
The dotted line represents the scalar field. 

The large logarithms are resummed by running v 
from ~ 1/r to ~ X/r in the solution of the renor- 
malization group equation. In our case, using that 
V$ = 1 at LL, we can solve the renormalization 
group equation for V s (for simplicity in the rest of 
this section we take the large N c limit): 



dv 



7T 



(18) 



with LL accuracy. We obtain (with this precision 
V s = E s ) 



E s = V S = -- 



(19) 



This result is correct with A ra+1 ln n A precision and 
is, obviously, gauge independent. 

We could consider the naive extrapolation to 
strong coupling of the above result. It does not agree 
with the AdS/CFT conjecture. We obtain A 2 ^ ver- 
sus the expected vA. This is not surprising as we 
have no reasons to believe our result should hold at 
strong coupling, the opposite limit with respect to 
which our result has been computed. 

We could consider playing a little bit and see how 
higher orders in the multipole expansion affect the 
result obtained in Eq. (Til))) . Therefore, we consider 
the following RG equation 1 



dv 



2±V S 

7T 



2 A 
— 7 
3tt 



• 2 K 3 ■ 



(20) 



1 Note that this result should only be considered for illustra- 
tive purposes, since contributions potentially of the same 
order (or even more important) have not been included in 
the anomalous dimensions. 



We obtain 

V. = 



A 1+2 ^ 



l + f(l-A^) 



(21) 



For this result is not possible to take the strong cou- 
pling limit, since it becomes imaginary. 

We would like to finish this section with some re- 
marks. The order at which infrared logarithms ap- 
pear in the static potential is easily visualized in the 
effective theory. In the case of QCD infrared loga- 
rithms first appear at 0(af/r). This follows from 
the multipole expansion suppression of the interac- 
tion of the octet and singlet field through the chro- 
moelectric field. There is no such suppression in the 
supersymmetric case, as we can see from the third 
line of Eq. ([5]), which explains the appearance of 
such effects already at 0(a 2 /r). On the other hand 
these effects can not come from gluons alone, since 
it would correspond to the pure QCD case. There- 
fore, the logarithms have to be associated with the 
scalars. This is seen quite clearly in the effective 
theory. 



IV. THE STATIC POTENTIAL AT ONE 
LOOP 

In this section we compute the singlet static po- 
tential and energy at one loop (as well as some 
other matching coefficients). This computation is 
necessary for an eventual complete NLL evaluation. 
Moreover, this example will allow us to visualize how 
the factorization between the soft and ultrasoft scale 
takes place in dimensional regularization including 
finite pieces. 

We first compute the soft piece, ie. the static po- 
tential. Diagrams with self energy cancel with di- 
agrams with internal vertices Therefore, only 
the ladder diagram (and the crossed one if we go 
beyond the large N c ) has to be considered. Note 
that this is only true in the Feynman gauge. In any 
case, the final result is gauge independent. There- 
fore, the expression of the bare static singlet poten- 
tial in D = 4 + 2e dimensions in momentum space 
reads 



Vs. 



or „ 2 1 Ji r 2,2. r[i-6]r 2 [ £ ] 

-2C F9 -jl-C^k 23+2t7r2+er[2e] 



-0(g 4 )} 



(22) 



After subtraction of the divergent piece in the MS 
scheme one can take the D — > 4 limit in the poten- 
tial. It reads 



V 



1 



A 



sMS = -2C F g 2 - {1-2- ln(k/u) + 0(X 2 ) 



(23) 
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We see that there is no finite piece associated to 
the logarithm in the MS scheme. Eq. (f2"3")) is the 
initial condition of the singlet static potential for the 
eventual renormalization group equation at NLL. At 
NLL we would also need the initial conditions for 
V$ and Va- Using arguments analogous to those in 
Ref. we can conclude that there are no one loop 
corrections to those matching coefficients and 



V* = 1 + C(A 2 ) . 
V* = 1 + 0(\ 2 ) 



(24) 
(25) 



V A = 1 + 0{\ 2 ) . (26) 
The bare potential in position space reads 
2C F g 2 r~ 2e ir[l/2 + e 



An 



r l/2+e 



1 



(27) 



2 3 7r 2 + e r[2e]r[l/2 + e] 



0(g 4 ) 



In the MS it reads 



2C F g 2 r- 2 * l T[l/2 + e} 
V sMs 7Z Z ZTFTTe l 1 ( 28 J 



An 



r n 



l/2+c 



2 3 7r 2 + e r[2e]r[l/2 + e] 



C A gW 1 



+ 7£ -ln(47r) +0(. 9 4 ) 



47r 2 \e 
which in four dimensions reduces to 



V s ,ms ~ 



2 -^l{l + 2±M V r) +lE ]+0{\ 2 ) 
r n 



Eq. (f2"9"| is the relevant object to be introduced in a 
Schroedinger equation. In order to obtain E s (r) we 
also need the ultrasoft contribution, ie. we need to 
compute the diagram in Fig. Q] (with the V$ vertex) 
including finite pieces. From this computation we 
obtain the following correction to the singlet static 
energy 



SV?J! = 2C F g 2 (V -V s ) 
+0(9'). 



1+2£ r[3 + 2e]r[-2-2e] 
2 2 + 2e 7r 3 / 2 + e r[3/2 + e] 
(30) 



Finally, the energy of two static sources in the fun- 
damental and anti-fundamental representation in 
J\f = 4 gluodynamics reads 

E s (r) = V S , B + 5V s u i = V Stm + SV s % , (31) 

which in four dimensions reduces to the following 
expression 

= _2C F a 1 f + 2 A +lE -i ]+ 0(\ 2 ) 

r y it 

(32) 



V. ORDINARY WILSON LOOP AND ITS 
ASSOCIATED STATIC ENERGY 

Besides the locally 1/2 BPS Wilson loop consid- 
ered in the previous sections, one could also consider 
the ordinary Wilson loop, which we define by elimi- 
nating the interaction of the static sources with the 
scalars in Eq (2), i.e. 



W c - T^TrPe 



-ig f c drA^x^ 



(33) 



One can then redo the analysis of the previous sec- 
tions for this case. Here we would like to compute 
the leading logarithmic corrections to the static en- 
ergy and potential. In order to do so we can use 
the effective Lagrangian in Eq. ijHJ) eliminating the 
terms proportional to the scalars fields 2 , rescaling 
by a factor 1/2 the static potentials in Eq. ([9]), and 
redefining A = N c a s /2. This implies that, to the 
order of interest, 71 = and 73 is 1/2 the value 
quoted in Eq. (|17p . From this exercise, we learn the 
important lesson that the leading In a s correction to 
the ordinary Wilson loop static energy in JV = 4 su- 
persymmetric Yang-Mills theory is the same to the 
one obtained in QCD [7|. The difference with QCD 
is due to the fact that a s itself does not run. This 
implies that there are no C(a™ +3 ln n (a s )) terms in 
the static singlet potential except for n = 1 (unlike 
in QCD [O]])- These findings can be summarized in 
the following equation, which is correct with NNLL 
accuracy, 



(29) E.(r) = - 



C F a s 



1 + aia s + a2<X 



(34) 

The coefficient ai = N c /ir has been obtained in Ref. 

We confirm this result. The coefficient a 2 is at 
present unknown. 

For the static singlet potential one should replace 
Ca& s by v jr in the logarithmic term, where v would 
correspond to the factorization scale. 



VI. CONCLUSIONS 

Eqs. fT9"|), |(52l and (03) are the main results of 
our paper. 

We have obtained the singlet static energy (and 
potential) with LL and NLO accuracy for the 1/2 



2 At present, we can not claim that those coefficients vanish 
at any order in a s but, at most, they are O(a^). This 
implies that the contribution from those terms is suppressed 
compared with those due to the Tr{OES} term. 
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BPS static Wilson loop. We have provided with ex- 
pressions at arbitrary dimensions, as well as with 
the formalism (based on effective field theories) , that 
could be relevant for future computations. We have 
illustrated how the merge of the soft and ultrasoft 
contribution takes place in the case of the singlet 
static energy including finite pieces. This may be of 
help in order to visualize how things will work in the 
case of QCD, where this mixing takes place at three 
loops. 

A naive interpolation to strong coupling of the 
renormalization group improved LL result does not 
agree with the supergravity conjecture. 

For the ordinary Wilson we have computed the 



singlet static energy and potential with NLO accu- 
racy and we have also performed the resummation 
of logarithms with NNLL accuracy up to the initial 
matching condition a.2- 
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